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Abstract. Working in the framework of a nonrelativistic quark model we evaluate the spectra and semilep- 
tonic decay widths for the ground state of doubly heavy S and Q baryons. We solve the three-body problem 
using a variational ansatz made possible by the constraints imposed by heavy quark spin symmetry. In 
order to check the dependence of our results on the inter-quark interaction we have used five different 
quark quark potentials which include Coulomb and hyperfine terms coming from one-gluon exchange, plus 
a confining term. Our results for the spectra are in good agreement with a previous calculation done us- 
ing a Faddeev approach. For the semileptonic decay our results for the total decay widths are in a good 
agreement with the ones obtained within a relativistic quark model in the quark-diquark approximation. 
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1 Introduction 

Even though only recently the mass of a baryon with two 
heavy quarks has been measured experimentally \, these 
systems have been being studied for more than a decade. 
Working with a system with two heavy quarks one can 
take advantage of the constraints imposed by heavy quark 
spin symmetry (HQSS). This symmetry amounts to the 
decoupling of the heavy quark spins in the infinity heavy 
quark mass limit. In that limit one can consider the total 
spin of the two heavy quarks subsystem to be well defined. 
This result, that we shall assume to be valid for the actual 
heavy quark masses, will simplify the solution of the three- 
body problem. 

In this contribution we shall present results for masses 
and total semileptonic decay widths. We have also an- 
alyzed other static observables as well as form factors, 
differential decay widths and angular asymmetries of the 
weak decays. For a detailed account of the full calculation 
see Ref. E] 

In Table \I\we summarize the quantum numbers of the 
baryons considered in this study. 



2 The model 

Once the centre of mass (CM) motion has been removed, 
the intrinsic Hamiltonian that describes the inner dynam- 
ics of the baryon is given by 



E 



Hf + V hlh2 (ri - r 2 , spin) - 



Vi • v 2 



Table 1. Quantum numbers of doubly heavy baryons analyzed 
in this study. S, J p are strangeness and the spin parity of the 
baryon, / is the isospin, and is the spin parity of the heavy 
degrees of freedom. I denotes a light u or d quark . 
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where r%, r 2 are the relative positions of the hi, hi heavy 
quarks with respect to the light quark q, M = m/^ + 



.l/m h . + l/m q 



and = d/d r , , j 



M 



m h2 + m q 

1,2. Vhjq and Vh 1 h 2 are the heavy-light and heavy- heavy 
interaction potentials. Note the presence of the Hughes- 
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Eckart term that results from the separation of the CM 
motion. 

For the quark quark interaction we have considered 
five different phenomenological potentials, one suggested 
by Bhaduri and collaborators |3] (BD) and four sug- 
gested by B. Silvestre-Brac and C. Semay HflS] (AL1, AL2, 
API y AP2). All of them include Coulomb and hyperfine 
terms coming from one-gluon exchange and a confining 
term, and differ in the form factor used for the hyperfine 
term, the use of a form factor in the one gluon exchange 
Coulomb term or in the power of the confinement term. All 
free parameters had been adjusted to reproduce the light 
and heavy-light meson spectra. Details on the potentials 
can be found in Refs. 

For the interactions considered, the total spin and in- 
ternal orbital angular momentum commute with the in- 
trinsic Hamiltonian, and thus are well defined. In this work 
we will study the ground state of baryons with total angu- 
lar momentum J = 1/2, 3/2 so we can assume the orbital 
angular momentum to be 0. This implies that the spatial 
wave function only depend on r\, r 2 and r\ 2 = |f i — r 2 \. 
We will also assume that taking the total spin of the heavy 
degrees of freedom to be well defined, as obtained in the 
infinite heavy quark mass limit, is a good approximation. 
That will allow us to write the wave function in a simple 
way (see Ref. for details). 

The spatial part of wave function will be determined 
using a variational method in which we will assume the 
following functional form: 

3C 2 (ri,r 2 ,r 12 ) = N F B (r 12 ) <j> hiq (n) <t> h2q (r 2 ) (2) 

where N is a normalization constant, 4>hjq is the S- Wave 
ground state wave function ip.j(r.j) of the single particle 
Hamiltonian H^ p corrected at large distances: 



(t>h 3q {r.j) = (1 + a s Tj) <Pj(rj), j = l,2 



(3) 



The heavy-heavy Jastrow correlation function F B will be 
given as a linear combination of gaussians: 



F B (r 12 ) 



4 

£< 



-6j(ri2+d 3 ) 2 
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where on, a; i ^ 1, bi and di are free variational parame- 
ters. The values that we get for the variational parameters 
are compiled in Ref. 0- 

We have also used the wave function obtained in this 
model to study different doubly 5(1/2+) -> B'(l/2+) 
baryon semileptonic decays involving a b — ► c transition 
at the quark level. We have worked in the spectator ap- 
proximation with only one-body currents. 
The differential decay width reads 



dr = 8\V cb \ 2 m B ,G 2 
x5 4 (p — p' — k 



d 3 P ' 
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d 3 k' 



F {2ir) 3 2E' B , {2n) 3 2E i?l (2tt) 3 2^ 
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(5) 



where \V c b\ is the modulus of the corresponding Cabibbo- 
Kobayashi-Maskawa matrix element, ra^ is the mass of 



Table 2. Doubly heavy £ masses in MeV. 
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the final baryon, Gf is the Fermi decay constant, p, p', 
k and k' are the four-momenta of the initial baryon, final 
baryon, final anti-neutrino and final lepton respectively, 
and C and TL are the lepton and hadron tensors. 
The lepton tensor is given as 

C^(k, k') = k^k a + k! a k^ - g" a k ■ k' + ie^ al} k' a kp(<o) 

where we use the convention e 0123 = —1, = (+, —,—,—) 
The hadron tensor is given as 



* (0)7mU- 75)^(0) 



B,r p 



B',r' p' 



V (0) 7ff (7-75)f 6 (0) fl,rp (7) 



with \B,r p) (\B' ,r' p 1 )) representing the initial (final) 
baryon with three-momentum p (p') and spin index r 
(r r ). The baryon states are normalized such that 

(r p | r 1 p') = (27r) 3 (£(p)/m) 5 rr , S 3 (p - p') (8) 

We compute the widths similarly as we did in Ref. |6] for 
baryons with a heavy quark. 



3 Results and discussion 

The mass of the baryon is simply given by the expecta- 
tion value of the intrinsic Hamiltonian. In tabled we give 
our results for doubly heavy E baryons, while in table 
are the results for the doubly heavy Q ones. Our central 
values correspond to the results obtained using the AL1 
potential, while the errors quoted take into account the 
variations found when using the other potentials. That 
also applies to the quoted results for Ref. 4 , obtained 
with the same interaction potentials but within a Faddeev 
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Table 3. Doubly heavy O masses in MeV. 
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Table 4. Semileptonic decay widths in units of 10 1 GeV. We 
have used \V c b\ = 0.0413. I stands for I = e, n 
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approach. When comparison with this work is possible we 
find an excellent agreement between the two calculations. 
Besides we give predictions for states not considered in 
the study of Ref . [4] . We also compare with other theoret- 
ical models. All calculations give similar results that vary 
within a few percent. From the experimental side only 
the mass of the S cc has been measured. The experimental 
value for M~ cc obtained by the SELEX Collaboration 
is 100 MeV smaller than our result. Note nevertheless that 
no account is given of the systematic error. There are also 
lattice calculations, by the UKQCD Collaboration [7], of 
the masses of the doubly charmed S cc , S* c , D cc and i?* e 
baryons. Our results are within errors of the lattice deter- 
minations. 

In table H we present our results for the semileptonic 
decay widths for the different processes under study. Our 



central values again correspond to the results obtained us- 
ing the AL1 potential, while the errors show the variations 
when using the other four potentials. The biggest varia- 
tions appear for the BD potential, with differences of the 
order of 7 ~ 12%. We compare our results with the predic- 
tions of different models. For that purpose we need to fix 
a value for \V c b\ for which we take \V c b\ = 0.0413. Our re- 
sults are in reasonable agreement with the ones in Ref. [Tl] 
where they use a relativistic quark model evaluated in the 
quark-diquark approximation. For r(£i, c — > S cc ) we also 
agree with the value of Ref. J7| obtained using heavy 
quark effective theory. A much smaller value for the same 
width is obtained in the relativistic three-quark model 
calculation of Ref. ^S]. In Ref. ^Sj, where they use the 
Bethe-Salpeter equation applied to a quark-diquark sys- 
tem, they obtain much larger results for all transitions. 
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